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Abstract 

The ground-state properties of two-dimensional liquid ^He at zero tem- 
perature are studied by means of a quadratic diffusion Monte Carlo method. 
As interatomic potential we use a revised version of the HFDHE2 Aziz poten- 
tial which is expected to give a better description of the interaction between 
helium atoms. The equation of state is determined with great accuracy over 
a wide range of densities in the liquid phase from the spinodal point up to 
the freezing density. The spinodal decomposition density is estimated and 
other properties of the liquid, such as radial distribution function, static form 
factor, momentum distribution and density dependence of the condensate 

fraction are all presented. 
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I. INTRODUCTION 



In recent years a great deal of interest has been devoted in the hterature to the study 
of quantum boson hquids in restricted geometries. From a theoretical point of view thin 
films of liquid ^He adsorbed on different solid substrates have been studied using variational 
techniques, based on the hypernetted-chain/Euler-Lagrange (HNC/EL) theory,^"^ Density 
Functional theories (DFT),^~^ and also Monte Carlo techniques. ^^'^^ According to these 
investigations, liquid ^He films display different behaviours depending on the strength of the 
substrate potential. 

In the case of weak-binding and long-range potentials, such as some alkali metal sub- 
strates, the helium film does not exhibit at all two-dimensional characteristics. By adding 
^He atoms to the substrate surface, three-dimensional clusters form and helium uniformly 
covers the surface only for coverages many atoms thick (prewetting transition). For some 
substrates a stable coverage is found only for infinitely thick films (wetting transition). 

More interesting, from the point of view of the role of dimensionality, is the case of sub- 
strates with deep and narrow potential wells for which the degree of freedom perpendicular 
to the surface is practically "frozen out" and low-coverage films are stabilized at the surface 
forming two-dimensional (2D) systems. To this category belong graphite, sohd H2 on glass 
and some alkah metal substrates like Li and Mg. On these substrates the growth of the first 
liquid layers is predicted to proceed via layering transitions:^'^'"^ by increasing the surface 
coverage single atomic monolayers develop and become stable one on the top of the other. 
This prediction has been confirmed experimentally for helium on graphite by heat capacity^ 
and third sound measurements.^ 

Variational calculations of ground-state and dynamic properties of single ^He monolayers 
at zero temperature reveal a striking 2D behaviour over a wide range of coverages. ^'^ Such 
calculations seem to indicate that single monolayers of liquid helium adsorbed on strong- 
binding substrates represent physical realizations of 2D homogeneous quantum liquids. 

A condition which can be relevant for theories aiming to describe ^He films is that 



2 



they possess the correct hmiting behaviour of homogeneous three-dimensional and two- 
dimensional systems. Homogeneous three dimensional (3D) ^He is a well studied system 
both experimentally and theoretically. In the case of homogeneous 2D *^He no experimental 
results are available and "exact" results coming from ab initio calculations can be extremely 
useful. 

From the point of view of Monte Carlo simulations, ^He in confined geometries has 
already been the object of several studies. Variational Monte Carlo techniques have 
been applied to the study of both inhomogeneous films on substrates^^ and more extensively 
to homogeneous 2D **He, where the simulation is easier. However, these methods rely 
heavily on the choice of the variational wavefunction used in the simulation. 

Prom the side of more exact Monte Carlo methods, thin ^He films on molecular hydro- 
gen surfaces at low temperature have been recently studied by Wagner and Ceperley^^ using 
Path Integral Monte Carlo (PIMC) techniques. Strictly two-dimensional ^He has also been 
the object of Monte Carlo calculations. ^^'^^ At T = Whitlock et al}^ have calculated, using 
Green's Function Monte Carlo (GPMC), the equation of state and other ground-state prop- 
erties of liquid and solid ^He in two dimensions, giving an estimate of the freezing and melt- 
ing densities. At finite temperature, PIMC techniques have been employed by Ceperley et 
al}* to investigate the superfluid transition, which in 2D belongs to the Kosterlitz-Thouless 
universality class. 

In the present work we present a Diffusion Monte Carlo (DMC) calculation of homoge- 
neous liquid ^He in two dimensions at zero temperature. Our main purpose is to provide an 
exhaustive calculation of the equation of state and other relevant properties of the ground 
state such as distribution function, structure form factor, momentum distribution and con- 
densate fraction for a wide range of densities within the liquid phase. In particular, we have 
performed a detailed analysis of the equation of state in the region of negative pressures 
in order to give an accurate estimation of the spinodal density of the 2D ^He hquid. The 
spinodal density is defined as the density below which the system becomes macroscopically 
unstable against density fiuctuations. Its determination is relevant for the physics of ^He 
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monolayers on strong-binding substrates because the spinodal density coincides with the 
coverage at which the uniform monolayer breaks down into 2D clusters in equilibrium with 
the vacuum, and it is thus the density at which the first layering transition occurs. ^'^ 

In our DMC simulation we have used as interatomic potential a revised version of the 
HFDHE2 Aziz potential^^ which was proposed by Aziz et al}'^ in 1987 and is known as 
HFD-B(HE) potential. This renewed interatomic potential was recently used by two of us^^ 
(J.B. and J.C) to study 3D bulk helium at zero temperature and it has been proven to give 
results for the density dependence of the pressure and the system compressibility which are 
closer to experimental data than the older Aziz potential. We are thus inclined to think 
that it gives also better results for two-dimensional ^He. 

The structure of the paper is as follows: in Sec. II we outline briefly the quadratic 
Diffusion Monte Carlo method used to solve the Schrodinger equation. In Sec. Ill we 
present results on both the equation of state in the liquid phase and on the distribution and 
structure functions. In Sec. IV are collected the results for the density dependence of the 
condensate fraction, which differ substantially from the previous analysis by Whitlock and 
coworkers, and of the momentum distribution for different values of the density. A brief 
discussion and conclusions are included in Sec. V. 



II. COMPUTATIONAL METHOD 

The aim of Diffusion Monte Carlo algorithms is to solve the time dependent Schrodinger 
equation of a system of N particles in imaginary time 

5^(R,t) 



dt 



(if-i?)^(R,t) , (1) 



where R = (ri, r^r) is the configuration vector for the positions of the N particles and t 
is measured in units of h. In Eq. (|I]) the Hamiltonian H has the usual form 

H=-^Wl + V{K), (2) 
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and E is a parameter representing an energy shift. Provided the initial wavefunction 
\E'(R, t = 0) has a nonzero overlap with the ground state $0(1^)) the solution of Eq. 
(|I|) gives exactly the ground state wavefunction in the asymptotic limit of large times: 
^(R,t ^ 00) = $o(R). 

When numerically integrating Eq. importance sampling techniques are employed to 
guide quickly the solution towards the ground state. To this aim the Schrodinger equation 
(H) is rewritten for the function 

/(R,t) = V^T(R)vl>(R,t) , (3) 
where iPt(R) is a time independent trial function. One gets 

-D V^/(R, t) + D Vr(F(R)/(R, t)) + {El{R) - E) /(R, t) = i /(R, t) . 



9/(R,t) 



dt 

(4) 

Eq. (^) has the form of a classical diffusion equation for the distribution function /(R, t). 
The first term describes the diffusive process with diffusion constant D = /2m. The 
second term contains the drift force 

F(R) = 2V't(R)"'Vr^/'t(R) , (5) 

which drives the system towards the region in configuration space where the trial function 
'?/'t(R) is relevant. The third term in Eq. (H) represents a branching term which depends 
on the local energy -El(R) = ?/'t(R)~"^-^^t(R)- 

By introducing the time dependent Green's function 

^(R', R, At) = (R' I e"^^* I R) (6) 

the solution of Eq. can be written formally as 

/(R', t + At) = j dR G(R', R, At)/(R, t) . (7) 

If the Green's function G{R! , R, At) is known for infinitesimal time steps At, the asymptotic 
solution for large times /(R, t 00) can be obtained by solving iteratively Eq. (|^). 



In the most usual implementations of DMC algorithms to solve the many body 
Schrodinger equation, (^(R', R, At) is approximated up to order At for small time steps. 
In this case, the obtained energy of the ground state depends linearly on the time step At 
and several calculations for different values of At are needed in order to extrapolate the 
correct result in the limit At — 0. 

Recently a quadratic DMC algorithm^'' has been proven to work efficiently in the descrip- 
tion of ^He droplets^^ and bulk 3D liquid ^He.^^ This method relies on an expansion of the 
Green's function G(R', R, At) up to order (At)^, which generates a quadratic dependence of 
the energy eigenvalue on At, permitting thus to avoid the extrapolation to the limit At ^ 
by working with a single time step. More details concerning this procedure can be found in 
Ref. 18. 

In the present work we have used in all the calculations At = 2 x lO^^r, where r is the 
appropriate time unit: r = ma'^ jlYi' (a = 2.556A). Time step analysis have been performed 
at the density p = 0.275cr~^, close to the liquid equilibrium density, and at p = 0.400o"~^, 
near the freezing density, and no changes have been observed in the results for the energy 
eigenvalue by reducing the time step At. 

The mixed estimator {ipT \ O \ $0) of a generic operator O, is the direct output of DMC 
algorithms. If the operator O commutes with the Hamiltonian of the system, the mixed 
estimator coincides with the pure expectation value on the ground state ($0 I O \ $o)- Ifj 
on the contrary, the operator O does not correspond to a conserved quantity, its expectation 
value on the ground state can be obtained by means of a linear extrapolation^^ 

($0161 $0) = 2(^T I 6 I $0) - (^T I O I i^r) . (8) 

The above method to get pure estimators is the most widely used in DMC simulations which 
must thus be supplemented by a variational Monte Carlo (VMC) calculation to determine the 
variational expectation value {ipT \ O \ ipr)- The linear extrapolation (^), which is obtained 
by writing the ground state wavefunction $0 as '^'o(R) = '?/'t(R) is correct to linear 

order in the functional variation 6ip. 

6 



To go beyond this approximation, removing the dependence on the trial wavefunction, 
several algorithms have been proposed to calculate pure estimators. ^^'^^ On this line an 
algorithm based on "forward walking" has been recently presented by two of us^^ (J. C. 
and J. B.) which is easy to incorporate in the original Monte Carlo algorithm and allows 
the calculation of pure expectation values of coordinate operators 0(R) with satisfactory 
stability and reliability The method applied to bulk 3D liquid helium^^ gives values for 
the particle-particle distribution function and static structure factor in very good agreement 
with experimental results and without any significant dependence on the function used 
as importance sampling. The results given in the present work for the potential energy, 
distribution and structure functions are calculated by using this method. 

Another important parameter in the calculation is the population of walkers n^, which 
represents the number of points Rj in configuration space at which the distribution function 
/(R, t) is sampled. In our calculations we have used a mean walker population of riy, — 450. 
At the densities p — 0.275(T~^ and p — 0.400(T~^ the walker population was increased to 
riyj = 900 and no appreciable change in the results was observed. 

Our simulation box contains 64 particles. At the 2D saturation density po — 0.284(T~^ 
this corresponds to a simulation box length of ~ 38 A, roughly 14 times larger than the 
mean interparticle distance at this density. Such box size is large enough to neglect finite 
volume effects. In fact, in 3D bulk hquid ^He finite size effects have been proven to be 
negligible already for a box length of 8 times larger than the mean interparticle distance. 

Finally, as importance sampling we have used a simple McMillan^^ two-body trial func- 
tion 

..(R)^nexp(-i(4)). 

We have taken for all the densities b — 1.205a". This value, which minimizes the energy at 
the equilibrium density, coincides with the one found by Whitlock et al}^ in their VMC cal- 
culation. For the highest density calculated, p = 0.420(T~^, we have also used as importance 
sampling the two-body function proposed by Reatto^^ 

7 



,,,K)^nexp(-i(4)"-r-(-(^)l). (^0) 

with L = 0.2, A = 2.0(7, A = 0.6a and b = 1.225a. Although this trial function gives at the 
density p = 0.420cr~^ a VMC energy lower than the McMillan function (P), no appreciable 
change in the results of the DMC simulation is found. 

III. RESULTS 

In this section, we present our numerical results for the energy and structure properties 
of the ground state. First we analyze the equation of state in the liquid phase, comparing 
our results with the ones presented in Ref. 11. Ground-state properties such as the radial 
distribution function and the static structure factor are also discussed in the first subsec- 
tion. The second subsection contains our results concerning the density dependence of the 
condensate fraction and the momentum distribution. 

A. Equation of state and structure properties 

In Ref. 11, the equation of state for both the liquid and the solid phase of 2D homoge- 
neous ^He has been investigated by means of VMC and GFMC techniques. The interatomic 
potential employed in this study is the two-body IIFDIIE2 Aziz potential. The estimated 
value for the liquid freezing density from the GFMC calculation is pi = 0.443(T~^. For the 
sake of consistency we have calculated the energy per particle using the Aziz potential at 
two densities in the liquid phase: at p = 0.275a~^, close to the equilibrium density, and at 
p = 0.400cr~^, near freezing. The values for the energy obtained from our DMC calculation 
are in good agreement with the GFMC results reported in Ref. 11 for the same two densities. 

Once the equivalence between our DMC algorithm and the GFMC algorithm employed 
in Ref. 11 has been tested by using the same interatomic potential, we have proceeded to the 
study of the equation of state in the liquid phase with the HFD-B(IIE) potential proposed 
by Aziz et al.^'^ (henceforth referred to as Aziz II). In Table I we present the results of our 

8 



DMC calculations for the total, potential and kinetic energy per particle for some of the 
densities calculated. As discussed in Sec. II the values of the potential energy per particle 
have been obtained by employing the algorithm for pure estimators. 

The Aziz II potential is slightly more attractive than the Aziz potential and the binding 
energies at the different densities are therefore somewhat larger with respect to the ones 
obtained in Ref. 11. For example, at p = 0.275(t^^ near to equilibrium the energies per 
particle are E/N = -0.8519 ± 0.0044 K and E/N = -0.8950 ± 0.0019 K for the Aziz and 
Aziz II potentials respectively. In Fig. 1 we show our DMC results for the equation of state 
in the liquid phase together with the GFMC results of Ref. 11. 

The equation of state of 2D liquid ^He is usually fitted by using a polynomial cubic 
function of the form 



where e = E/N and po is the equilibrium density. In Table II we report the values of the 
parameters which best fit our results and we compare them with the values reported in 
Ref. 11. The values for the equilibrium density are very close, whereas appreciably different 
are the values for the B and C parameters. These differences affect the predictions for the 
density dependence of the surface pressure and compressibility, as well as the estimation of 
the spinodal density. The two fits are shown in Fig. 1 together with the Monte Carlo data. 
The cubic polynomial fit (|lT]) fits our data rather well and no significant improvement in 
the quality of the fit is found by increasing the order of the polynomial function used. 

Once the equation of state function e(p) is known, one can calculate straightforwardly 
the surface pressure, defined as -P(p) = p^{de/dp), and the isothermal compressibility 




(11) 




(12) 



From the inverse compressibility 1/k one obtains the velocity of sound 




(13) 
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In Fig. 2 we show the comparison between the surface pressure obtained from our equation 
of state and the result of Ref. 11. An appreciable difference is found in the regimes of low 
and high densities. A larger difference is found by comparing the predictions of the two fits 
for the velocity of sound, shown in Fig. 3. The velocity of sound at the 2D equilibrium 
density is c(po) = (92.8 ± 0.6) m/sec, which is nearly 3 times smaller than the velocity of 
sound at the saturation density of 3D bulk liquid ^He, c"^^(po) = 238.3 m/sec. 

At T = 0, and in the hmit of a 2D film, the third sound velocity C3 coincides with the 
velocity of sound c^^ of the purely 2D hquid.^ Then, for low coverages c^^ must reproduce 
quite accurately the C3 in monolayers adsorbed on strongly binding substrates. A direct 
comparison between c^^ and experimental C3 for low-coverage films is however difficult by 
the present uncertainties in the determination of the liquid coverages. 

The density dependence of the velocity of sound is linear over a wide range of densities. 
Only when approaching the spinodal density, where the system becomes unstable against 
density fluctuations, the velocity of sound drops suddenly to zero and the compressibility 
diverges. The estimation of the spinodal density from our fit gives: Psp = (0.228 ±0.002)(T~^, 
which is significantly smaller than the value Psp — 0.247(T~^ obtained from the data of Ref. 
11. The spinodal density of 2D homogeneous hquid ^He has been estimated also by Clements 
et al} from the fit of their HNC/EL variational results for the equation of state: they find 
Psp = 0.202(7-2. 

At the spinodal density Psp — 0.2280"-^ the mean interparticle distance is Vgp — 
{1/npspY^'^ — 3.O2A. For 3D bulk liquid ^He an accurate estimation of the spinodal density 
has been given in Ref. 30 with the result p^J^ — (0.264 ± 0.002)(7-^, corresponding to a 
mean interparticle distance of r^^ — (S/ATrplpY^^ — 2.47A. It is interesting to notice that 
at the 2D freezing density pi = 0.4430"-^, estimated in Ref. 11, the mean interparticle dis- 
tance is only r; = 2.17A, whereas at the 2D equilibrium density the interparticle spacing is 
Teq — 2.71A. The corresponding values of the average interatomic distance in the 3D case are 
respectively: rf^ = 2.10Aand r^^ = 2.22A. The mean interparticle distances of the 2D and 
3D liquid ^He barely overlap; at freezing the particles are only 3% further apart in 2D than 
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in 3D, at the equilibrium density the difference increases to 22% and the same difference 
persists down to the spinodal density. 

Important information on the structure of the ground state is obtained from the two-body 
radial distribution function 



9Kr^2)- /|$o(ri,...,r^)Pdri...dr^ ^ > 

and from its Fourier transform, the static structure factor 

5W^I.,/*e-teW-l)4<55|^. (15) 



with 

N 



p^ = Y^e^^-^. (16) 



i=l 



Both these quantities can be calculated using the method for pure estimators described in 
Ref. 25. In Fig. 4 we show the results obtained for the radial distribution function at 
three different densities. As the density increases more peaks at large interparticle distances 
appear. At the highest density shown p = 0.420cr~^, which is just before freezing, four peaks 
are clearly visible in the distribution function: a clear indication that the system is close 
to sohdification. In Table III, we report the position Vm and the height g{rm) of the first 
peak in the radial distribution function for some of the densities calculated. The height of 
the first peak in g{r) increases with the density and it shifts towards smaller interparticle 
distances. It is interesting to compare the height of the first peak in the radial distribution 
function for the 2D and 3D systems: in 2D at the equilibrium density g{r„i) ~ 1.25 whereas 
the corresponding value in 3D is 1.38. This is a clear indication that the 2D system is 
more dilute and possesses less correlations at equilibrium than its 3D counterpart. Close to 
the freezing density the heights of the first peak in g{r) of the 2D and 3D systems become 
comparable. In fact, as previously discussed, the mean interparticle distance of the two 
systems become similar when freezing is approached. 

In Fig. 5 we show the static structure factor for three values of the density. As the 
density increases the peak in S{k) increases and the values at the lowest momenta accessible 
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in our calculation increase. Due to phonon excitations the form factor S{k) is expected to 
go to zero in the long wavelength limit as S{k) ~ k/2mc. As the density decreases and 
the spinodal density is approached the velocity of sound c drops to zero and consequently 
the slope in S{k) diverges. This behaviour, which has been observed in the variational 
calculations of Ref. 1, agrees qualitatively with our ab initio calculations. 



B. Condensate fraction and momentum distribution 

Another quantity of great interest is the one-body density matrix p(r', r), which is related 
to the change in the ground-state wavefunction when a particle is removed from position r 
and replaced at position r'. For a homogeneous system p(r) is defined as 

where ■0(0) and ■0^(r) are, respectively, the field operators which destroy a particle from 
position r = and create one at position r. 

In the DMC algorithm the mixed matrix element ('0t|'0^(?^)'0(O)|$o)) involving the trial 
wavefunction ipT, can be calculated by averaging over the asymptotic distribution function 
/(R, t — > oo) the relative change in the trial wavefunction when a particle is displaced from 
position Fj to Fj -I- r 

{i^Tmr)m\<^o) = / /(R^ t PC) (^T(ri, . + r, r^)/0.(r„ r^)) dr,...dv. ^^^^ 

J /(f*-)^ ~^ oo)ari...arjv 

The calculation of the ground-state expectation value ($o|'*/'^(^)'*/'(0)|$o) can not be obtained 

straightforwardly from our method for pure estimators because it involves the knowledge of 

$o(R'')^o(R'), where R' is the configuration vector with one particle displaced by r, instead 

of <I>q(R). To calculate the one-body density matrix p(r) we employ thus the extrapolation 

technique described in Sec. II. For the highest density calculated, p — 0.420cr~^, we have 

used as importance sampling both the McMillan wavefunction (9) and the one proposed by 

Reatto (10). The extrapolation technique gives, within statistical errors, the same result for 

p(r) in the two cases. 
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The asymptotic limit of p(r) gives the fraction of particles uq condensed into the zero- 
momentum state 

rin = lim p(r). (19) 

r— »oo 

In Fig. 6 we show the results for the condensate fraction obtained at different densities 
ranging from the spinodal point up to the freezing density. We have fitted our data with 
the quadratic polynomial 

no(p) = no(po) + a + b ' , (20) 

where po is the equilibrium density po — 0.284(7"^. The values of the parameters giving the 
best fit are the following 

no(po) = 0.233 ±0.001 

a = -0.583 ±0.006 (21) 
6 = 0.44 ±0.02. 

The value for the condensate fraction at equilibrium density no(po) is consistent with the 
estimation reported in Ref. 14 no(po) — 0.22, obtained by extrapolating to zero temperature 
the PIMC results for the algebraic decay of the one-body density matrix. Our results for the 
condensate fraction are somewhat smaller than the ones reported in Ref. 11; for example 
at low density, p — 0.275(7"^, we find Uq — 0.251 ± 0.005 whereas the GFMC calculation 
of Ref. 11 gives n^^^^ = 0.36 ± 0.05. This discrepancy between DMC and GFMC results 
persists over the whole density range and only at very high density, p = 0.400(T~^, Hq'^'-^ 
and Uq^^^ become consistent. The reason for this is unclear, certainly it is not due to the 
revised version of the Aziz potential used in our DMC simulation. We have repeated our 
DMC calculation of the condensate fraction using the Aziz potential for the two densities 
p — 0.275cr~^ and p — 0.400cr~^ and no difference was found with the results obtained with 
Aziz II. It is interesting to notice that near the freezing density, where the mean interparticle 
distance in the 2D and 3D systems are comparable, also the condensate fraction is nearly 
the same (no 4%). 
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By Fourier transforming the one-body density matrix p(r), one gets the momentum 
distribution of the system 



In Fig. 7 we show the momentum distribution plotted as kn{k) for three different densities. 
As the density increases, more and more particles leave the condensate state and the max- 
imum in kn{k) shifts towards higher momenta. The shoulder present at high momenta is 
interpreted in 3D bulk ^He as coming from the zero-point motion of the rotons.^^ A similar 
interpretation could be valid also in the 2D case, in this case though the shoulder appears 
more pronounced than in 3D, particularly at low densities, and shifted towards smaller /c's. 



The properties of homogeneous 2D liquid "^Hc at zero temperature have been investigated 
by means of a quadratic diffusion Monte Carlo method. As interatomic potential we have 
used a revised version of the Aziz potential which is expected to be more accurate in de- 
scribing the interaction between helium atoms. The energy per particle has been calculated 
for a wide range of densities from the spinodal point up to the freezing density, providing 
an accurate determination of the equation of state function. The spinodal density of the 
2D system is estimated. The radial distribution function and the static structure factor are 
calculated for various densities employing a recently devised method for estimating pure 
expectation values. The fraction of particles in the condensate state has been calculated 
for various densities and the density dependence of the condensate fraction has been also 
estimated. 
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FIGURES 

FIG. 1. Equation of state for 2D liquid ^He. The solid circles correspond to the DMC energies 
obtained with the Aziz II potential (the statistical error bars are smaller than the size of the 
symbols); the solid line is the polynomial fit (11) to the calculated energies. The open diamonds 
are the GFMC results of Ref. 11 with the Aziz potential and the dashed line is the reported fit to 
these values. 

FIG. 2. Density dependence of the surface pressure as obtained from our fit to the equation of 
state function (solid line), and from the fit of Ref. 11 (dashed line). 

FIG. 3. Density dependence of the velocity of sound as obtained from our fit to the equation 
of state function (solid line), and from the fit of Ref. 11 (dashed line). 

FIG. 4. Radial distribution function for three densities: p = 0.275(7"^ (solid line), p = 0.320(7~^ 
(short dashed line), p = 0.420cr~^ (long dashed line). 

FIG. 5. Static structure function for three densities: p = 0.275cr~^ (solid line), p = 0.320cr~^ 
(short dashed line), p = 0.420(T~^ (long dashed line). 

FIG. 6. Density dependence of the condensate fraction. Solid circles with error bars: results of 
DMC calculations; solid line: fit from eq. (20). 

FIG. 7. Momentum distribution for three densities: p = 0.275(T~^ (soUd line), p = 0.320(7'"^ 
(short dashed line), p = 0.420(T~^ (long dashed line). 
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TABLES 

TABLE I. Results for the total and partial energies from DMC calculations. The potential 

energies per particle have been obtained from the calculation of pure expectation values. 



p(a-2) E/N{K) V/N{K) T/N{K) 

0.235 -0.8480± 0.0016 -3.850± 0.014 3.002± 0.014 

0.255 -0.8799± 0.0020 -4.234± 0.022 3.354± 0.022 

0.275 -0.8950± 0.0019 -4.680± 0.015 3.785± 0.016 

0.320 -0.8599± 0.0017 -5.722± 0.013 4.862± 0.013 

0.340 -0.7932± 0.0031 -6.213± 0.026 5.419± 0.026 

0.380 -0.5697± 0.0047 -7.390± 0.026 6.820± 0.026 

0.420 -0.1524± 0.0047 -8.532± 0.024 8.379± 0.025 



TABLE IL Parameters of the equation of state (11) for our DMC results and the GFMC results 

ofRef. 11 

Parameter DMC GFMC 

/9o((T-2) 0.28380 ± 0.00015 0.28458 

eo{K) -0.89706 ± 0.00061 -0.8357 

B{K) 2.065 ± 0.014 1.659 

C{K) 2.430 ± 0.035 3.493 

X^/iy 0.99 1.45 
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TABLE III. Position Vm and height g{rm) of the first peak in the radial distribution function 

(A) g{r^) 

0.235 4.09 1.209 ± 0.005 

0.255 4.03 1.233 ± 0.005 

0.275 3.98 1.253 ± 0.002 

0.320 3.87 1.323 ± 0.003 

0.340 3.84 1.363 ± 0.009 

0.380 3.63 1.426 ± 0.013 

0.420 3.61 1.522 ± 0.006 
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